Dynamics and delocalisation transition for an interface driven by a uniform shear flow 
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We study the effect of a uniform shear flow on an interface separating the two broken-symmetry 
ordered phases of a two-dimensional system with nonconserved scalar order parameter. The inter- 
face, initially flat and perpendicular to the flow, is distorted by the shear flow. We show that there 
is a critical shear rate, 7 C cx l/L 2 , (where L is the system width perpendicular to the flow) below 
which the interface can sustain the shear. In this regime the countermotion of the interface under 
its curvature balances the shear flow, and the stretched interface stabilizes into a time-independent 
shape whose form we determine analytically. For 7 > ~/ c the interface acquires a non-zero velocity, 
whose profile is shown to reach a time-independent limit which we determine exactly. The analytical 
results are checked by numerical integration of the equations of motion. 



There is growing interest in understanding the effects 
of shear on the dynamical properties of statistical systems 
far from equilibrium The most extensively studied 
case is perhaps the approach to equilibrium of a system 
quenched below its critical point: domain growth is heav- 
ily influenced by the external shear and new dynamical 
exponents appear The determination of these expo- 
nents, together with the problem of the validity of dy- 
namical scaling, are the most challenging tasks in this 
context ||. 

The transition from a disordered, high-temperature 
phase to an ordered, low-temperature one, however, is 
not the only context where an applied shear may play a 
major role. An interesting alternative problem is to in- 
vestigate cases where the shear may by itself introduce 
a novel dynamics in a state otherwise ordered and sta- 
ble. In the case of spinodal decomposition of a binary 
fluid, for example, we may let the system evolve until 
a stable, entirely separated state is reached, and then 
apply a shear flow normal to the interface between the 
two phases. Further evolution will then occur, with a 
competition between the shear, which tries to stretch the 
interface, and diffusion of the two constituents, which 
tends to straighten it. 

More generally, a natural question in this context is to 
what extent the applied shear is able to perturb the sta- 
ble initial state. More specifically, we may ask: Is there a 
critical value for the intensity of the shear, beyond which 
the system is unable to restore itself in a stable stationary 
state ? 

In this paper we shed some light on the above questions 
by studying the dynamics of a flat interface subjected to a 
transverse shear flow. The motion of such an interface in 
the case of conserved dynamics has been recently studied 
in 0| , where the existence of a critical value of the shear 
beyond which stationarity was lost was not reported: for 
the shear rates studied, the stretched interface always 
reached a stable steady state. In the present work we 
study the deterministic dynamics of a similar interface 
under shear, but for the case of nonconserved dynamics, 
in a system described by a scalar order parameter. Ising- 



like systems, such as twisted nematic liquid crystals [gj, 
display a behavior that can be described by this model. 
For weak shear we find, similarly to Q| , that the interface 
reaches a stationary profile, in which the curvature forces 
acting on the interface compensate the shear. In this way 
the interface slips relative to the moving boundaries, and 
acquires a steady-state profile with no net velocity in the 
laboratory frame. However, for shear rates larger than 
a critical value we find that the interface cannot sustain 
the strain and a time-independent state can no longer be 
reached. In this regime, the interface departs indefinitely 
from its initial condition and becomes delocalised. In- 
stead of a stationary spatial profile, the interface acquires 
a stationary velocity profile. It should be noted that the 
speed of the interface at the boundaries is always smaller 
than that of the boundaries themselves, i.e. the contact 
points still slip with respect to the moving boundaries, 
but not enough to keep the interface stationary in the 
laboratory frame. 

The system we study is a two-dimensional strip, 
bounded in the y direction between the values — y and 
y , so its width is L = 2y Q . We consider a uniform 
shear flow, given by the simple velocity profile v = 7?/e x , 
where 7 is the shear rate and e x is a unit vector in the 
x-direction, which is thus the direction of the flow. At 
time t = 0, the interface is given by the equation x = 0, 
i.e. it is a flat segment connecting the boundaries and 
separating two regions with opposite value of the order 
parameter. We will consider here only the zero tempera- 
ture dynamics of such a system. 

As a physical boundary condition we assume that the 
interface is perpendicular to the boundaries at the points 
of contact: any different condition would create an infi- 
nite force at such points due to the uneven curvature. 
We can understand better this argument in the context 
of the Ising model, by considering an interface which is 
not normal to the boundary: the contact spin on the 
acute angle side has an excess of neighbouring spins with 
opposite sign, and it therefore flips. This makes the in- 
terface become normal to the interface. This argument 
only holds if the microscopic flipping time tq is smaller 



1 



than the time needed by the flow to shift the spin, that 
is if 7 < 1/tq. We will find a critical value of the shear 
7c -C 1/tq and therefore we can consistently assume an 
interface normal to the boundaries at the contact points. 
Moreover, we have checked that an Ising simulation with 
free boundary conditions produces an interface which is 
indeed normal to the boundaries at the contact points. 

The deterministic dynamics of a non conserved order 
parameter, </>(x, £), is described by the time-dependent 
Ginzburg-Landau equation 0], 



(i) 



where r Q is the relaxation time for order parameter fluc- 
tuations in the bulk, £ G is the interfacial width, and V(<p) 
is a symmetric double well potential. Starting from Eq. 
(|l|) one can show using standard methods that the nor- 
mal velocity of an interface separating the two phases is 
proportional to the local curvature. This is the Allen- 
Cahn equation R|, v = —D V • n, where V • n is the 
curvature, D = Co/ T o is the diffusion constant for fluctu- 
ations in the bulk phases, and n is the local normal to 
the interface. The Allen-Cahn equation is actually more 
general than the specific dynamical equation ([!]) , since it 
simply expresses the dynamics of an interface driven by 
its surface tension. 

In the absence of shear, a flat interface is stable. When 
a shear is applied, however, the Allen-Cahn equation has 
to be modified: 



v = —D V • n + 72/ e„ 



(2) 



where the additional term represents the advection of 
the interface by the flow, leading to a distortion of its 
initially planar form. Describing the interface profile by 
the function x(y,t), which gives the displacement in the 
flow direction as a function of the position in the direc- 
tion transverse to the flow, it is easy to show that the 
divergence of the normal is given by: 



V ■ n 



d 2 y x 



while the velocity, v, of the interface can be written as 

d t x 



(3) 



[l + (d v xf]^ ' 



(4) 



where dtx is the velocity of the interface in the direction 
parallel to the flow. Combining (||), (||) and (^) we obtain 
an equation for the interface profile x(y,t): 



df.x = D 



1 + (dyX) 



(5) 



Prior to further analysis, it is convenient to rescale space 
and time, 



X = x/y Q , Y = y/y a 



tD/yl 



(6) 



where now Ye [—1,1]. Note that in this way we are also 
introducing a rescaled interface velocity V — vy a /D. In 
terms of the rescaled variables, Eq. (||) reads 



d T X 



d\X 



1 + (d Y X) 



aY 



IVI/D 



(7) 



Due to the rescaling, all the dimensionfull parameters 
have been absorbed into the dimensionless effective shear 
rate a. Given the geometry of the problem, we clearly 
expect X(Y, r) to be an odd function of Y. We will there- 
fore limit our analysis to the domain Y € [0, 1]. 

In order to find a stationary solution X S (Y) of equation 
(Q), we set d T X = 0, to obtain (where primes indicate 
derivatives with respect to Y) 



X" 



1 + {X' S ? 
with boundary conditions 



-aY 



X' S (1) = = X S (0), 



(9) 



where the first condition follows from the requirement 
that the interface be perpendicular to the boundaries, 
and the second from the fact that X(Y) is an odd func- 
tion. Integrating once, and imposing the boundary con- 
dition at Y = 1, gives 



X' S (Y) = tan 



-(1-F 2 



(10) 



A second integration, incorporating the boundary condi- 
tion at Y = 0, gives the stationary interface profile, 



X{Y) 



tan 



r(i-^ 2 ) 



(11) 



This interface profile is plotted in Figure 1 for different 
values of the effective shear rate. 

From equation ( |To| ) we see that the function X S (Y) has 
a maximum in its derivative at Y = 0. When a = -k the 
derivative diverges at this point indicating that the inter- 
face is parallel to the system boundaries. A value a > tt 
would imply a discontinuous derivative X' S (Y) and thus 
an interface profile X S (Y) with two cusps: such a profile 
would be unphysical, because of the infinite curvature at 
the cusps. We conclude that for a > tt no stationary 
solution is possible. This implies that there is a critical 
value of the dimensionless shear rate, 



i.e. a critical shear rate 



4ttD 



(12) 



(13) 
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FIG. 1. Stationary interface profile as given by equation 
(p"l|), for three different values of the effective shear rate a. 
The initial configuration of the interface is given by the flat 
X = segment. 



X s (l), of the asymptotic time-independent solution as a 
function of the shear rate, in order to compare our an- 
alytic result, Eq. ([ll]), with the numerical integration of 
Eq. (0). The agreement is good, and it can be seen that 
it improves as the continuum limit is approached. The 
conclusion is that the stationary state given by Eq. ( pT| ) 
is indeed an attractor for the interface dynamics. 
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For a > a c (or 7 > j c ) we must find a time-dependent 
solution of the full equation (Q). In terms of the micro- 
scopic parameters of the systems, we have 



7c 



4< 2 1 



TQ 



(14) 



and therefore j c <C 1/tq as long as the interfacial width 
is much smaller than the size of the system, a condition 
that obviously holds for any reasonable system. There- 
fore, as anticipated in the Introduction, our assumption 
of an interface normal to the boundaries at the contact 
points is fully justified in the totality of the interesting 7 
regime. 

For a — > a c the contact point X(l) goes to infinity 
and the interface is no longer localized in a finite region 
of space. We have, 



X(l) 



a 



(15) 



In order to obtain this result it must be noted that for 
a c the integral in equation (11) is completely con- 



o 

centrated near the origin, z — 0, and it therefore de- 
pends very weakly on the value of Y , as long as Y is 
non- vanishing. This means that for a — > a c not only the 
contact point X(l), but all the points of the interface 
X(Y) with non-zero Y are found at the same position 
given by @. 

Before considering the a > a c regime, we want to check 
that the time-independent solution we have found is in- 
deed an attractor for the interface dynamics, when start- 
ing from the initial configuration X(Y,0) — 0. To do 
this, Eq. (0) was discretized and integrated numerically. 
In Figure 2 we plot the position of the contact point, 



12 3 4 

a 

FIG. 2. Numerical integration of equation (Q) (symbols) 
and analytical result (full line). The numerical integration 
approaches the analytical result as the spatial discretization 
scale, Ay, decreases. 

When a > a c the time-independent equation for the 
interface has no physical solution and the interface must 
move with a nonzero velocity parallel to the boundaries. 
A reasonable assumption is that, in the large-time limit, 
the velocity of every point along the interface is time- 
independent, i.e. d T X(Y,r) = V{Y,t) -> Voo(Y). Such 
a velocity profile for the interface implies that the non- 
linear term in equation (|?]) must be time-independent in 
the asymptotic limit. This condition must be satisfied in 
two different ways, according to the different regions of 
the y-domain. 

• Y ~ : From the time-independent equation (|^) 
we see that X' s '(0) = for all a < a c . It is natural then 
to assume that even in the phase a > a c the interface 
is asymptotically flat close to the centre, such that the 
nonlinear term vanishes in this region and the solution of 
the equation is trivial, 



X(Y,t) = aYr 



(16) 



This solution clearly satisfies Eq. (^), and the boundary 
condition X(0,t) = 0. However, it cannot be correct 
close to the boundary at Y = 1 because it does not rat- 
ify the boundary condition 0yX(1,t) = 0. We will find 
that it is correct in a domain < Y < Y\ , with Y\ to be 
determined below. 
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• Y ~ 1 : At Y = 1, the interface must satisfy 
the boundary condition cV_X"(l,r) = 0, and therefore 
we need to keep the nonlinear term in the equation of 
motion. However, as we already noted, this term must 
become independent of time for large times. This can be 
achieved by requiring that 



X(Y,r) = f(Y) + V T 



(17) 



so that the portion of the interface close to the system 
boundary moves with a velocity V a independent of Y. 
This solution also satisfies Eq. (0). We shall find that it 
is correct in the domain Y\ < Y < 1. Integrating equa- 
tion (0) with d T X(Y,r) = V a and dyX(l,T) = 0, we 
obtain 



dyX = tan 



~(i-y 2 )-v (i-y) 



(18) 



The two solutions eqs. ( |16| ) and ( |18| ) must be matched 
at the border between the two different regions. Let us 
define region I: < Y < Y\ and region II: Y\ < Y < 1. 
The matching point Y\ will be calculated later. We can 
write the complete interface velocity profile, for r — > oo, 



V X {Y) 



uY 

V = aYx 



I 

II 



(19) 



For the Y derivative of the interface profile, for r — > oo, 
we have 



tan[f(i-y 2 )-y (i-y)] 



I 

II 



(20) 



Furthermore, note that the value of V is known, from 
Eq. (|l9|), once Yi is determined. In order to calculate Yj 
we have to match the solutions in the different regions 
at this point. Note that the derivative, dyX, has to be 
continuous at this point to avoid an infinite curvature. 
From ( p0| ) , we see that this requires the argument of the 
tangent in region II to be equal to ir/2 when Y = Y\, 
This gives 



Y x = 1 



Vo = a 1 - 




(21) 



Note that the velocity V of the interface at the contact 
point, Y = 1, is smaller than the velocity of the flow at 
this same point which, in our dimensionless variables, is 
equal to a. Our analytical solution was tested by nu- 
merical integration of the dynamical equation, and the 
results are shown in Figure 3. Numerical and analytical 
results are fully consistent. 

The convergence of the interface velocity at the con- 
tact point, d T X{l, t) — > V for t — * oo, was found to be 



very slow. The function d T X(l,t) is monotonically de- 
creasing in time: initially it is equal to the flow velocity 
at the contact point a, but eventually it must decrease, 
in order to develop a curvature to keep the interface per- 
pendicular to the boundary. A closer inspection of Figure 
3 shows that for any finite time the cusp in the velocity 
profile is rounded out and this function is smooth. 
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FIG. 3. The velocity profile of the interface as a function 
of Y, with a — 12. It can be clearly seen that for large times 
the velocity profile converges towards the function Voo(Y) of 
eq. ©. 

However, for very large times, it is not unreasonable 
to approximate the actual velocity profile with the same 
piecewise linear form as in eq.(po|) (see Figure 4). Of 
course, we have to introduce a time-dependent value, 
Yi(r), for the matching point, with Yi(r) — > Y\ for 
t — ► oo . Within this approximation we can therefore 
write 



d Y X(Y, T ) 



tan[f(l-Y 2 )-y o (r)(l-y)] 



/ 

II 

(22) 



where V {t) = aYi(r) and Yi(r) = Yi + SYi (r) . The new 
value of the matching point Y\{t) is fixed, as usual, by 
imposing the continuity of 3yX(Y, t). This gives 

tan-V") « ZI - — „2[i_y l(T )]2 (23) 
I ar I 

Setting Yi(t) = Y X + 5Y\ and recalling that f (1 - Y{) 2 = 
^ from the definition of Yi in equation ( pl| ) , we get 



5Y X 



a 2 (l-Yi)r ' 



(24) 



and for V(1,t) 



4 



V(1,t)^V + 



1 1 



(25) 
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FIG. 4. The nature of the approximation leading to equa- 
tion f0). 



In this paper we have analytically and numerically 
studied the zero temperature dynamics of an interface 
subject to a transverse shear flow, in the case of non- 
conserved dynamics. We find a critical value, a c = n, 
of the dimensionless shear rate, a = jy^/D, beyond 
which a steady state cannot be reached and the interface 
moves with a constant velocity. In terms of the thickness, 
L = 2y a , of the sample, this gives 7 C = AirD/L 2 , i.e. the 
transition occurs at lower shear rates for wider systems. 

While no such critical value was reported for the con- 
served case analyzed in [Q , it is unclear to us why there 
should be a major difference in this respect between con- 
served and nonconserved dynamics. Indeed, in the lat- 
ter case, the dynamics is slower, due to the conservation 
constraint, and the system should be less capable of sus- 
taining the shear than in the nonconserved case. We 
conjecture, therefore, that a similar transition occurs for 
systems with conserved dynamics. 



In order to check this last result numerically, it is con- 
venient to compute the time derivative of V(l,r), which 
does not contain the constant V a . We have 



d T V(l,r) 



1 



1 



(26) 



This last result is tested in Figure 5, where we may see 
that again numerical integration of the dynamical equa- 
tion and analytical results coincide. 
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FIG. 5. The quantity d T V(l,r) as a function of r in a 
log-log plot. The full line has slope —2, for a comparison 
with the analytical result eq.(^). 
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